Polaronic effects in electron shuttling 



m 
o 

(N 
C 



G. A. Skorobagatko/'B I. V. Krive, 1 ' 2 and R. I. Shekhter 2 
1 B. Verkvn Institute for Low Temperature Physics and Engineering of the National 
Academy of Sciences of Ukraine, 47 Lenin Avenue, Kharkov 61103, Ukraine 
Department of Physics, University of Gothenburg, SE-412 96 Goteborg, Sweden 

(Dated: January 31, 2013) 

Shuttle-like mechanism of electron transport through a single level vibrating quantum dot is con- 
sidered in the regime of strong electromechanical coupling. It is shown that the increment of shuttle 
instability is a nonmonotonic function of the driving voltage. The interplay of two oppositively 
acting effects - vibron-assisted electron tunneling and polaronic blockade - results in oscillations of 
the increment on the energy scale of vibron energy. 
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I. INTRODUCTION 

The modern trends in miniaturization of electronic de- 
vices eventually led to fabrication of single molecular 
junctions and molecular transistors (see e.g. review^). 
The electric properties of single molecular transistors 
(SMTs) in many cases are similar to the analogous char- 
acteristics of single electron transistors (SETs) fabricated 
in two-dimensional electron gas. SMTs demonstrate such 
effects as Coulomb blockade and Coulomb blockade os- 
cillations on gate voltage. The significant difference be- 
tween SMT and semiconducting SET is that the former 
can function even at room temperatures that makes them 
to be very promising basic elements for nanoelectronics. 

Another specific feature of molecular transistors is the 
interaction between electronic and vibrational degrees 
of freedom. The electron in the process of tunneling 
through the molecule can excite (and absorb at finite 
temperatures) molecular vibrational quanta (vibrons) - 
the phenomenon known as " phonon-assisted tunneling" 2 . 
The opening of inelastic channels results in appearance 
of additional peaks (side-band peaks) in differential con- 
ductance. For weak electron-vibron interaction the mag- 
nitudes of inelastic peaks are much smaller then the value 
of the elastic resonance peak. The situation is changed 
in the regime of strong electron-vibron interaction when 
nonperturbative (and, in particular, polaronic) effects de- 
termine electron transport through a vibrating molecule 
(see e.gii). 

Polaronic effects are most pronounced in the case when 
the molecule (quantum dot) is well-separated from the 
leads and the width, T, of conducting molecular states is 
small compared to other energy scales (temperature T, 
driving voltage eV). In this case the mechanism of elec- 
tron transport through the vibrating molecule is (inelas- 
tic) sequential tunneling via the real polaronic interme- 
diate state. The amplitude of this tunneling is exponen- 
tially suppressed since the wave functions of free electron 
in the leads and polaronic state (Holstein polaron) in the 
dot are almost orthogonal. This effect (named as Frank- 
Condon 3 , 4 or polaronio^ blockade) strongly suppresses 
elastic channel of electron transport and changes the 



temperature behavior of conductance. Recently, Frank- 
Condon blockade was observed in experiment of elec- 
tron tunneling through a suspended single-wall carbon 
nanotube^. 

A one more novel phenomenon appears for vibrating 
quantum dot when the matrix element of electron tun- 
neling to the left and to the right lead differently de- 
pends on the position of the dot center of mass. This 
is always the case when the dot (molecule) vibrates in 
the direction of electron tunneling and then the effect 
of electron shuttling takes place at finite voltages^. In 
papers-,— the problem of electron shuttling was consid- 
ered for a model of single level (sq) vibrating quantum dot 
weakly coupled to the leads of noninteracting electrons. 
It was shown that in the regime of weak electromechani- 
cal coupling the shuttle instability occurs at bias voltages 
eV > 2(eo + fo^a) (ftwo is the vibron energy) and the in- 
crement of instability is tq ~ (To/H)X\t- Here To is the 
level width, A is the dimensionless electron-vibron inter- 
action constant and A* = xo/h (xq is the amplitude of 
zero-point fluctuations of quantum dot, It is the electron 
tunneling length). Both coupling constants A, At were 
assumed to be small in Refs^^. 

In the problem of electron shuttling the electron-vibron 
interaction strength X(V) linearly depends on the driv- 
ing voltage. Therefore, at sufficiently high voltages the 
regime of strong electron-vibron interaction (A > 1) is 
realized. In this regime polaronic effects could play sig- 
nificant role in electron shuttling. In the present paper 
we reconsider the problem of shuttle instability for the 
regime of strong electromechanical coupling. 

We derived the equation of motion for the shuttle av- 
erage coordinate x(t) assuming only the weak character 
of dot-lead interaction. It was shown that the increment 
r s (V) of shuttle instability is a nonmonotonic function 
of bias voltage with a maximum at eV m ~ huio{d/ 'xq) 2 , 
where d ^> xq is the distance between the leads. The 
maximum value of the increment is sensitive to value of 
coupling constant At. The interplay of two effects - the 
increase of r s caused by the increase of inelastic chan- 
nels which contribute to the increment when rising the 
applied voltage and the decrease of r s with the increase 
of bias voltage caused by polaronic blockade - results in 
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FIG. 1: Schematical diagram of Single Electron Transis- 
tor (SET) with vibrating quantum dot (QD). Here A(eV) 
denotes the "shifted" (due to polaronic shift) bias voltage- 
dependent fermionic level; huo is the energy of vibrational 
mode; Tl(r)(x) are the coordinate-dependent tunneling am- 
plitudes and ±eV/2 are the chemical potentials of the leads 
ef = 0. Characteristic distances in the QD: d is the distance 
("gap") between the leads; l t is the tunneling length of the 
electron in the QD; x(t) is average coordinate of the shuttle. 



oscillation of r s on small energy scale fia>o- 

Our results show that polaronic effects determine the 
physics of electron shuttling in the case of moderate 
or strong mechanical damping when the transition to a 
shuttle-like regime of electron transport is possible only 
at sufficiently high bias voltages. 



II. THE MODEL 

Our starting point is the model of vibrating quantum 
dot weakly coupled to the leads of noninteracting elec- 
trons. This model was repeatedly considered in the liter- 
ature for the problem of electron transport in molecular 
transistors (see e.gi and referencies therein). We ex- 
pand this model to the problem of electron shuttling^ by 
taking into account the explicit dependence of tunneling 
amplitude on the center of mass coordinate of quantum 
dot. 

For simplicity we will study the case of a single level 
(with the energy £o) quantum dot coupled to a single 
vibronic mode (with the energy ftw^). The total Hamil- 
tonian of our system is 
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j=L,R 



(./) 



where 



H 



k 



Ek 



(1) 



(2) 



is the standard Hamiltonian of noninteracting electrons 
(efe) in the left (j = L) and right (j = R) leads, fij is 
the corresponding chemical potential: — ^r — eV 
( V is the driving voltage); (a/cj) are the cre- 

ation(destruction) operators. The Hamiltonian of vibrat- 
ing quantum dot takes the form (see e.gJ£) 



Eqc + c — Ei(b + + b)c + c + huj b + b . 



(3) 



where Ei is the characteristic energy of electron- vibron in- 
teraction (see below), c + (c) and b + (b) are fermionic and 
bosonic creation (destruction) operators with commuta- 
tion relations {c, c + } = 1; [6,6+] = 1. The tunneling 
Hamiltonian is 



\ij(X c .m.)a>kj( 



h.c. 



(4) 



In Hamiltonian Eq.([J| we take into account the depen- 
dence of tunneling amplitude tj on the center of mass 
coordinate X c .m. of quantum dot. In quantum descrip- 
tion the coordinate X c/m , becomes an operator 



X, 



xo 



V2 



P 



(5) 



where xq = ^h/Muo, (M is the mass of QD) is the 
amplitude of zero-point oscillations of quantum dot. In 
Eq.([5]) we defined also the dimensionless momentum 
operator P with the canonical commutation relation 
[X, P] — i. In what follows the tunneling amplitude is 
model 8 - by the exponential function 



t 3 {X) = t oj exp(jX t X), j = (L, R)=(-, 



(6) 



Here At = xo/lt and It is the tunneling length. 

Notice that electron- vibron interaction term in Eq.Q 
in our model originates from the electrostatic interaction 
of charge density on the dot with the electrostatic poten- 
tial produced by the leads 8 . It is convinient to charac- 
terize this interaction by the dimensionless bias voltage- 
dependent coupling constant 



A = 



eV 

2hLUn 



^ « 1 (7) 



(V is the bias voltage, d is the distance between the 
leads) . We use the notations usually accepted in the liter- 
ature on molecular transistors. Notice that our notations 
for the coupling constants differ from Refs<2r— . 

The problem of electron shuttling in the model Eqs.(fT])- 
§S§ was studied in Refs^^ for the case of weak electrome- 
chanical coupling (in our notations: A < 1, At < 1). 
In molecular transistors the electron- vibron interaction 
can be strong A > 1 (see e.gpi). Here we reconsider the 
problem of shuttle instability^ in the regime of strong 
coupling. 

To study electron transport in the presence of pola- 
ronic effects (A > 1) it is convenient to use unitary trans- 
formation (see e.gJ^-) which eliminates electron-vibron 
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interaction term in the dot Hamiltonian Eq.(|3]). Shut- 
tle instability results in appearance of classical time- 
dependent coordinate x(t) of quantum dot. It means 
that bosonic operators b and b + acquire c-number part 
a(t): 



b = a{t)+b; b+ =a*(t)+b+ 



(8) 



here b + (b) are the bosonic creation(destruction) op- 
erators, which describe quantized vibron modes 

(b) = (b + ) = ((. . .) denotes the thermal average). 

We transform total Hamiltonian (1) using the unitary 
operator (see e.g.—) 

U = cxp(iXpn) , h = c + c , p = —j= (i + — bj . (9) 

The unitary transformed Hamiltonian takes the form 

H = UHU- 1 = e(t)c + c+haj b + b+ ^ H ( t 3) + h[ ]) 

j=L,R j=L,R 

(10) 

where e(t) = e - [Xx(t) + \ 2 /2}huj , x(t) = V2Rea(t). 
The time- independent shift (~ A 2 ) of the dot level 
is called polaronic shifti The transformed tunneling 
Hamiltonian in Eq. dTUl) is 



(i) 



E 



a+e- iX Pc+h.c. 



*0iE [ij{X t )a+V^c + h.c. 



(11) 



Here X t is the transformed coordinate operator 
X t = UXU- 1 = X + Xh and 

Vj = exp(jXtX + iXp) , 
toj = too exp(j^) , j = (L, R) = (-, +) . (12) 



III. EQUATIONS OF MOTION 



Here we denote by Hj the tunneling Hamiltonian defined 
by Eq. fTTj) . In our model equations ([Hf take the form 



dX 
~dt 



= u P-X E Jj(X,P) 



3=L,R 



^ = ~^ X ~i E 3Hj(X,p), (16) 

j=+/R,-/L 

where Jl,r are the current operators 
. i j 



J, 



E[ 6 ^-^ 



(17) 



These operators satisfy (as it should be) the continuity 
equation 



dn 
~dt 



Jl + Jr. ■ 



(18) 



With the help of Eq. (TT8")) we can rewrite the first expres- 
sion of Eq. (|16[) in the following form 



P = LO 



and the second equation in (1161) is transformed to 



(19) 



- 2 {x + XcM) + .lX 



E 3Sj(X,p), 



j=-/L,+/R 



(20) 

where p = P — p(t). 

To make the operator equations (Eq.(fT6t or Eqs. (fT9|) ), 
(1201) complete we have to write down the equations of 
motion for fermionic operators afcj(a^ ) and c(c + ). These 
equations 



~~dT 



-i{e k - Hj)a kj (t) - it 0j V+{t)c(t) , (21) 



The Heisenberg equations of motion for dimensionless 
operators of coordinate X and momentum P 



X = -^=[h' +b) =.,■(/) 



P 



V2 



b+ - b) =p(t)+p. 



(13) 



(where x(t) — y/2Rea(t), p(t) = \/2Ima(t) ) can be 
represented in the form of Hamilton equations 

dX _ 1 dH dP _ 1 dH 

with the Hamiltonian H given by the following expression 
H=^(X 2 + P 2 )+ £ H,(X,p). (15) 

j=L,R 



dc 
dt 



-ie(t)c(t)-i E *oj^jf(*) a fcj(*) 



(22) 



k,j=L.R. 



are linear and they can be readily solved (the equa- 
tions of motion for creation operators are obtained 
from Eqs. (f2~Tj) , (|22l) simply by taking Hermitian conju- 
gation). We will follow Refsi^ and find the solution of 
Eqs. (|2"T]) . (l2"2"j) in the so called "wide band approxima- 
tion" (" WBA")(see e.gJ^L). The only difference of our sys- 
tem of equations (j2"Tj) . (|22l) from the corresponding one in 
Ref4 is the presence in our case bosonic operator factors 
Vj and Vy 1- . Formally, these factors make the level width 

f (see the definition below in Eq. (|24|) . which appears in 
Eq. (|2"2")1 after substituting in this equation the solution of 
Eq. (l2~Tj) . to be an operator as well. The dynamical equa- 
tion for the dot level operator c in WBA takes a simple 
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form 



takes the form 



dc 
dt 



e(t) 



At) 



-i(s k -Hj)t 



k,j=L,R 



where 



-f(f)= J2 



(23) 



(24) 



j=L,R 



Here Tqj — 2irpj\toj\ 2 and pj is the density of states 
in the j = L, R electrode (pj is assumed to be energy 
independent in the wide band approximation). 

To proceed further we have to make additional sim- 
plifications. Our purpose in this section is to derive 
the equation of motion for the average coordinate in 
the presence of fluctuations. Notice that in the Hamil- 
tonian Eq. (|10p electron- vibron interaction appears only 
in the tunneling term which for a weak tunneling can 
be treated perturbatively. In perturbation theory on 
the bare level widths Toj vibrons are decoupled from 
fermions and when evaluating averages of the products 
of fermion (F) and boson (B) operators we can use the 
decoupling procedure (FB) (F) g + jj t (B) g b . The av- 
erages of fermion operators are taken with "fcrmionic" 
part of the total transformed Hamiltonian (including the 
tunneling part) and averages of boson operators are cal- 
culated with the quadratic Hamiltonian of nonointeract- 

ing vibrons Hb = hujob + b. This decoupling procedure is 
the starting point of calculations in many papers dealing 
with the polaronic effects in electron transport in molec- 
ular junctions (see e.g. review^ and referencies therein). 

The linear equation (|2"3")l can be formally solved (us- 
ing time-ordering procedure) for the operator level width 
T(t). However the closed equation for the average coordi- 
nate x(t) can be obtained only in perturbation theory on 
Toj. In Ref.- such equation was obtained in the limit of 
weak electromechanical coupling A <C 1, At <C 1. Here we 
derive the equation for classical coordinate x(t) valid also 
for the regime of strong electron- vibron coupling A > 1. 

In perturbation theory on the bare level width Tqj we 
can neglect the time dynamics of the level width and 
replace T(t) by constant Tq — Tql + Tor. It is con- 
vinient in what follows to represent "vertex" operator Vj 
as a product of classical Tj{x(t)} and quantum Qj parts 
V 3=T] {t)Q 3 



o 



c(t) 

t r ft 

dt\Tj{t)Qj exp <i / dti 



k,j=L,R 



£k - Hj - e(*2) + 



(26) 



The equation of motion for the classical coordinate x(t) 
in perturbation theory on T^j can be readily obtained 
from the exact operator equation (|20[) by taking the av- 
erage and using the discussed above " fermion-boson" fac- 
torization procedure 



dt 2 



[x(t) + XN{x(t)}}+cu 2 x(t) 



^oA t 



E 3Hj{x{t)} 
(27) 



where 



N{x(t)} = (c+(t)c(t)) , Hj {x(t)} = (Hj(X,p)) 



It is useful to rewrite Eq. (|27p in terms of a new vari- 
able x v (t) — x(t) + XN{x(t)}. Notice that according to 
Eqs. (fTTj) . (l2l))) both quantities in Eq. (l2"Tl) - averaged tun- 
neling Hamiltonian Hj and the level occupation num- 
ber N{x(t)} - are proportional to the bare level width 
Hj,N{x(t)} oc T j. Since Eq.(J27|) is derived in the Born 
approximation (up to the second order in the tunneling 
amplitude) we can replace x(t) by x v (t) in the averaged 
quantities. Then Eq.([2~T|) for the variable x v (t) (shifted 
coordinate) takes the form of the Newton's equation de- 
rived in Ref£ 



dt 2 



where 



x v (t) + uj 2 x v (t) = F q {x v {t)} , (28) 



Fq {x(t)} = Xco 2 Q N{x(t)}-^ £ jH 3 {x(t)} (29) 



■ -/L 



This equation in the limit A <C 1, A t <C 1, when we 
can omit operator factors (Qj) in the vertex function 
Vj , exactly coincides with the corresponding equation for 
the coordinate of "classical" shuttle (see Refs.— ^). One 
can analyze shuttle instability by using either Eq. (|27[) or 
Eq. (l28|) . The only difference is the starting (equilibrium) 
position of the shuttle. This is x — for Eq. (|2"T|) and 
x v = XN ^ for Eq.® (see Ref£). We will use Eq.®. 



, . r x -/ m A . n , , » , ^ IV - THE INCREMENT OF SHUTTLE 

Tj(t) = cxp[jX t x(t)\ , Qj(x,p) = cxp(jX t x + iXp) , INSTABILITY 

j = (L,R) = (-,+) 

(25) The conditions for a shuttle instability can be found by 

analyzing linearized equation of motion. We will follow 
In the discussed approximation the solution of Eq. (|23[) Ref^ where these conditions were obtained for the regime 
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of weak electromechanical coupling. The linear integral- 
differential equation for the shuttle coordinate x(t) in the 
dimensionless units (t — >• uj$t, energy scale is Hloq) reads 
(see Appendix) 

,2 f +°° 

y{t) + m=- e -i E f m 

tfe/j (e) 1 2/m [X y (e, t)] AG^ +) (e) - jA 

(30) 

In Eq. (l3T)l) the following notations are introduced: 
f\y = r j exp(jAAj) 



Xij(s,t) = / dra;(T 
Jo 



■ exp 



and 



(e - fij) — (A — l)± iF /2 



(31) 



(32) 



where A = eq — A 2 /2 is the polaronic shift. Vibrational 
degrees of freedom result in Franck-Condon factors 



F l3 (fi) = exp[-(A 2 -A 2 )(l + 2n b )] 



X + jXt 



I A - j\t 

■h(2\\ 2 - \ 2 t Wn b {l + n b ))e~M 2 (33) 

Here h{x) is the modihed Bessel function of the sec- 
ond kind and n b — l/(exp/3 —1), (j3 = huJo/ksT), 
fj(e) = l/[exp(/3(e— fj,j))+l] arc Bose-Eistein and Fermi- 
Dirac distribution functions. Remind that all energies in 
above expressions are dimensionless (in the units of huo). 
In the limit A t = the vibron-induced correlation factor 
Eq. (l33p coincides with the well-known in the literature 
expression (see e.g. Refs. ( 10 ' 12 ). 

The solution of Eq. (l30l) with the initial condition 
x(0) = is 



x(t) = A)e rat sin(t). 
where Aq is an arbitrary constant and 



(34) 



= E 

+ /R 

E - 

■ -/L 



f E F«C8)(A-iAt) a /i(A-l + l) 

/ — — OO 



- Oj 



+ 00 

E 

1 — — OC 



Fy(/3)(A + iA t ) 2 /,(A-/-l) 



(35) 



is the dimensionless increment (when r s > 0) of shuttle 
instability. Our purpose here is to find conditions (driv- 
ing voltage) for the realization of shuttle motion in the 
presence of strong electron- vibron interaction. 



For a weak electromechanical coupling A C 1 , A f < 1 
we can neglect quantum and thermodynamical fluctua- 
tions of vibrons (they are of higher orders on coupling 
constants) and omit all terms in the right-hand side of 
Eq.® but I = 0. In this limit A ~ e , f 0j ~ T oj , 
Fi oj ~ 1 and Eq. (|35[) is transformed to the corresponding 
equation of Refs/'-. The increment Eq. ([33)l in this case 
coincides with the one found earlier Ref£ 



-(0) 



= E ^-j^fjieo + l) 



•_-/£ 



E ^(A + jA f ) 2 /,(eo-l) (36) 



(there is a misprint in Eq.(24) of Rcf. 9 - the sign-changing 
factor j is missing in the coefficients in front of distribu- 
tion functions). 

It is qualitatively clear that shuttle instability is most 
pronounced at low temperatures (3 ^> 1 when there are no 
thermally activated vibrons. At low temperatures shut- 
tle instability results in a sharp (step-like) features in the 
I—V curve of single electron transistor—. Finite tempera- 
ture effects broaden the transition region and atT> Hu}q 
the shuttle instability is less pronounced (nevertheless as 
far as r s > the instability results in a strong increase 
of electric current). In what follows we will consider only 
low-T effects and (for simplicity) the case of a symmetric 
junction (T 0L = T 0R = T , fi L = -fx R = eV/2, we set 
Fermi energy of the leads Ef = 0). At low temperatures 
the / > terms in Eq. (|35[) are exponentially suppressed 
in comparison with the negative I (formally, due to the 
factor exp(— 1/3/2) since I-i(x) — h{x) for the integer I). 
Physically, positive I corresponds to vibron absorption - 
an energetically forbidden process at T <C fiio . Nega- 
tive / describes emission of vibrons and the summation 
over I at finite bias voltage is limited by a certain l m (see 
below) . 

By using the well-known asymptotics of the Bessel 
function at small arguments Ii(x — > 0) <~ (z/2) l /l\ and 
replacing Fermi distribution functions in Eq. (|35p by the 
Heaviside theta-functions we obtain the desired formula 
for the increment (r s > 0) of shuttle instability 



.. .. =i4^exp(-A 2 + A 2 -AA t ) E 

1=0 



2 

rxx t 



(A + A t ) 
2! 



2/ 



= 2 „ _ 1} , cx P (AA t + 2A 2 )r (l m , (A + A t ) 2 ) , (37) 

where T(a, x) is the incomplete gamma function (see 
e.gfi^) and 



eV ( X\ 



(38) 



([x] denotes the integer part of x). We find 
from Eqs. (l3"Tj) . (|38[) that in the regime of weak 
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electromechanical coupling the instability occurs at 
eV > eV c ~ 2(e + an d the increment 

r s (A < l,A t < 1) ~ r = rAA t /2 (see Ref£) is a lin- 
ear function of bias voltage. Remind that we consider 
the case when the uncertainty in the quantum dot ini- 
tial position due to quantum fluctuations (~ xo) is small 
compared to the geometrical size d of the junction (that 
is our parameter r ( i = 2xo/d <C 1, Eq.([7])l. In this case at 
the threshold voltage V c the electron-vibron coupling is 
small X(V C ) = r d (e + 1)^1 (we always can put e = 
at the resonance condition) and polaronic effects arc not 
pronounced. Nevertheless there is a small negative cor- 
rection to the threshold voltage due to polaronic shift 
V c = 2(eo — A 2 (V c )/2 + 1) and the multiplicative renor- 
malization of the bare level width by quantum fluctua- 
tions r -> r exp(-A 2 + A 2 - AA t ) ~ T(l - A 2 - AA t + A 2 ) , 
as one can see from Eqs. (|37|)(l38)l . 

The increment r s is described by Eq. ([57|) at bias volt- 
ages in the interval eV c /2 < eV/2 < 2r7 2 . At higher 
voltages the electron distribution function of the right 
electrode (biased by — eV/2) in Eg. ([55)) Jr = 1 and 
the processes of inelastic electron tunneling to the right 
bank start to contribute to r s . Their contributions ac- 
cording to Eq. (l35l at low temperatures are negative 
and they could only diminish the increment. We show 
now that due to polaronic (Franck-Condon) blockade 
both the "right" and "left" contributions at voltages 
eV > 2r~ j 2 1 are exponentially suppressed and shuttle 
instability takes place in the finite interval of bias volages 
(we restore here the dimensions) 

2(e + Hu) ) <eV < 4tuj /rj (39) 

The finite series on I in Eq. (l37[) can be approximated 
as follows 

With the help of asymptotics Eq. ([40|) we ob- 
tain the following formula for r s at eV <C r^ 2 
(the corresponding electron-vibron coupling "constant" 
A = r d (eV/2) « rj 1 ) 

FA A 

r s ~^exp(AA t + 2A 2 ) (41) 

For a very high biases eV 3> r7 2 it is easy to show 
from Eqs. (f3"T)) . (l4"(Jl) (using Sterling's formula to estimate 
asymptotics of l m \) that r s oc exp[— A 2 (I — ln2)/2] — > 
at A 3> 1. So we see that the dependence of increment 
r s on the bias voltage is strongly nonmonotonic with the 
maximum at eV ~ r^" 2 . 

It is interesting to notice here that at the excitation 
energy Ed ~ huia(rd)~ 2 (the corresponding number of ex- 
cited vibrons Id — r^ 2 3> 1) the characteristic width, w, 
of the wave function of harmonic oscillator which repre- 
sents quantum dot in our model (Eq.©) is of the order of 
gap, d, between the leads (w ~ xq\/ (x 2 ) ~ XQ^/Td ~ d). 
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FIG. 2: "Weak" shuttle instability. The increment of shuttle 
instability (in the units of r/ftuo) as a function of bias volt- 
age for A t = xo/h = 0.2; Td = xo/d = 0.2 (solid line) and 
= ksT/Hujo — 0.2. The dotted line represents the result 
of Re^ extended to the region of strong electromechanical 
coupling. 
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FIG. 3: "Strong" shuttle instability. All parameters are the 
same as in Fig. 2, but for the coupling constant At = 0.7. 



The approximation Eq.(|40p does not reveal the fine 
structure (on the Hloq scale) of the dependence r s (V) . On 
this small energy scale one could expect the appearance 
of steps each time the additional vibration channel con- 
tributes to the increment (see Eq.([3"5|)). At low voltages 
the steps are slightly modified by the Franck-Condon fac- 
tors. With the increase of voltage in the regime of strong 
coupling A > 1 each additional channel modifies r s by 
the value of the order (few times smaller) of r s (sharp 
big steps). However between two sequential steps r s is 
diminished (due to the exp(— X 2 (V)) factor in Eq.([57P) 
approximately by the same amount. That is on the scale 
of vibron energy there are oscillations of r s . These oscil- 
lations are smeared out at temperatures T 3> HajQ. 

The dependence r s (V) is plotted in Figs. 2, 3 for the dif- 
ferent values of parameters of our model (for numerical 
calculations we used the basic Eq. (|3"T|) and put the inverse 
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temperature (3 — 10). For both figures we set r<j = 0.2, 
since the further decrease of r c i yields the results which 
practically coincide with Ref.— in the considering region 
of applied voltages. The dimensionless increment r s of 
the shuttle instability is measured in the units of dimen- 
sionless r = T/huo (eV is in the units of hu)o). In Fig. 2 
we set A t = r<2 = 0.2, that corresponds to the case when 
2l t ~ d. In Fig. 3 we put A t = 0.7 > r<j, that corresponds 
to the case when 2l t < d. The straight line on both fig- 
ures shows the increment for the "classical" shuttle of 
Refi^ for the same values of all parameters. 

We see from the figures, that, the dependence of in- 
crement on bias voltage is strongly nonmonotonic. The 
increment r s grows until eV < eV m ~ r^ 2 and decreases 
when V > V m . The magnitudes of r s in Fig. 3 (At = 0.7) 
are much greater than the ones in Fig. 2 (At = 0.2). This 
means, that shuttle instability is very sensitive to the 
value of tunneling lengths At. So, when l t <C d (r^ <C 1), 
the greater is At the stronger is shuttle instability (i.e. 
it develops on a shorter time scale ~ fj 1 )- In the 
presence of strong mechanical friction, characterized by 
phenomenological damping term, ijfx(t), introduced in 
Ea.(|2"7). and for 2l t < d (i.e. r d < 1), when A t < 1, the 
increment at all voltages could be less then the friction 
coefficient r s < 7/. In this case shuttle regime of electron 
transport is not realized. 



V. SUMMARY 

In this paper we have studied the influence of quantum 
and thermodynamical fluctuations on shuttle instability. 
These fluctuations are significant in the case of strong 
electromechanical coupling that can be realized in elec- 
tron transport through vibrating quantum dot at high 
bias voltages. 

It was shown that the increment of shuttle instabil- 
ity is a nonmonotonic function of bias voltage V with a 
maximum at eV m ~ hu>o(d/xo) 2 , which corresponds to 
the region of strong electron- vibron coupling \{V m ) ^> 1. 
At higher voltages polaronic blockade suppresses shuttle 
instability. The maximum value of increment is sensi- 
tive to the electromechanical dimensionless coupling con- 
stant At = Xo/h (h is the electron tunneling length). In 
the presence of mechanical friction (characterized by phe- 
nomenological friction coefficient 7/) and when At -C 1 
the increment r s ( V) at all voltages could be less then fric- 
tion r s (V) < 7/ and shuttle regime of electron transport 
is not realized. In the regime of strong electromechan- 
ical coupling At ~ 1, A 3> 1 the pumping of energy in 
mechanical shuttle motion could overcome even strong 
damping. 

We showed that in this regime the increment of shut- 
tle instability strongly oscillates on the energy scale of 
vibron energy Hloq. If the friction coefficient is compa- 
rable with the amplitude of increment oscillations, the 
small change of bias voltage (AV < fouo/e) drastically 
changes the regime of electron transport through a vi- 



brating quantum dot (from a phonon-assisted tunneling 
to a shuttle regime of electron transfer and vice-a- versa). 
One can speculate that reentrant transitions to a shuttle- 
like regime of electron transport will result in unusual 
current- voltage characteristics with pronounced negative 
differential conductance. 
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VII. THE APPENDIX 

To study the shuttle instability we linearize the Eq. (|28|) 
(Xx(t) <C 1, Xtx(t) <C 1). In the linear approximation 
we have Tj(x(t)) = t 0j (l + jX t x(t)), j = -/L;+/R and 
c(t) = co(i) + ci(i, x(t)), where 



(*) = ~i E toja k j(0)e 



-i(s k -Hj)t 



dt\Qj{t\) exp < i 



k,j=L,R 



e k — fij — A + -T 



(t - h) 



(42) 



and 



c x (t) = -i X t oj a k j{0)e 

k,j=-/L,+/R 

ft 

dti 



—i(ek—Hj)t 







j\ t x(ti) + A / dt 2 x(t 2 ) 
ti 



■Qj(h) exp < i 



e k - [ij - A + -r 



(t-ti) 



(43) 



Here A = £0 — A 2 /2 is the polaronic shift. 

Then, the linearized force F q (t) in the right-hand side 
of the Eq. (|28|) is represented as the sum of two contribu- 
tions: F q = F + Fx, where 

F = X(c+(t)c (t)) 

k,j=-/L,+/R 

E 

k,j=-/L. + /R 



jXttojic+^akj^e-^-^Qjit)) 



(44) 
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and 

F 1 =-X t x(t) Y, h 3 (Qj(t)a+(0y^-H)tc o{t)) 

k,j=L,R 

-X t x(t) Y M#(*)M°) e_<(e,, ~ w)1 £i(*)} 

k,j=L,R 
k,j=-/L,+/R 

£ iA t «o i (c+(i)a fej (0)e- i ^-^)*4(i)) 

k,j=-/L,+/R 

+\[(c+(t)c (t)) + (c+(t)c 1 (t))} 

(45) 

To proceed further, we calculate the averages in 
Eqs. (|4"4"l) . (|4"5j) . In perturbation theory on the level width 
r the averages of boson and fermion operators factorize 

= (a^m^miQjih)^)) , (46) 

where (o^.(0)a fc / 3 -(0)) = fj(e k )S kk ,. For non- 
interacting equilibrated electrons in the leads 
fj{ £ k) — [exp(,$(£fc — fj,j)) + is the Fermi-Dirac 
distribution function. 



It is easy to calculate boson correlation function in 
Eq. (|46j) . assuming vibrons to be at equilibrium at tem- 
perature T (this is a plausible assumption for a weak 
tunneling regime we are dealing with). The standard 
calculations (see Ref^) results in 

(Qt(ti)Qj(t 2 )) = £ FuWexpQUti-tn)), (47) 

l — — oo 

where Fij(f3) is defined in the main text (see Eq. ([33|) ). 

Now, by substituting Eq.([47]l into the Eqs. (j44|) . f|45]) 
and taking all time integrals in the limit t To, we 
obtain the right-hand side of Eg. ([25]) (i.c.F g ) as a linear 
functional of x(t). In this case the "force" term Fq does 
not depend on time 

F a = A(c+(t)c (t)) = \ E E FiMfA^-l) 

j— L,i? /— — oo 

= const 

(48) 

It determines the initial position of QD. The linear in x(t) 
term F\ takes the form of the right-hand side of Eq. (|30| . 
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